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ABSTRACT 

We use the Fourth Data Release of the Sloan Digital Sky Survey to investigate the relation 
between galaxy rest frame u — r colour, morphology, as described by the concentration and 
Sersic indices, and environmental density, for a sample of 79,553 galaxies at z < 0.1. We split 
the samples according to density and luminosity and recover the expected bimodal distribution 
in the colour-morphology plane, shown especially clearly by this subsampling. 

We quantify the bimodality by a sum of two Gaussians on the colour and morphology 
axes and show that, for the red/early-type population both colour and morphology do not 
change significantly as a function of density. For the blue/late-type population, with increasing 
density the colour becomes redder but the morphology again does not change significantly. 
Both populations become monotonically redder and of earlier type with increasing luminosity. 
There is no significant qualitative difference between the behaviour of the two morphological 
measures. 

Motivated by their long-standing use in astronomy and their ability to utilise information 
not necessarily used by the concentration and Sersic indices, we supplement the morpholog- 
ical sample with 13,655 galaxies assigned Hubble types by an artificial neural network. We 
find, however, that the resulting distribution is less well described by two Gaussians. There- 
fore, there are either more than two significant morphological populations, physical processes 
not seen in colour space, or the Hubble type, particularly the different subtypes of spirals 
Sa-Sd, has an irreducible fuzziness when related to environmental density. 

For each of the three measures of morphology, on removing the density relation due to 
it, we recover a strong residual relation in colour. However, on similarly removing the colour- 
density relation there is no evidence for a residual relation due to morphology. Therefore, 
either the morphology is not directly affected by the environmental density beyond the corre- 
lation to colour, or a single galaxy 'type' does not capture sufficient information. 

Key words: cosmology: observations - methods: data analysis - methods: statistical - galax- 
ies: fundamental parameters - galaxies: statistics 



1 INTRODUCTION 

The connection between the morphology of a galaxy and t he den- 
sity of its environment (e.g. IOemleilll974l : lDresslei][T98Qh . mea- 
sured in a variety of ways, is an important clue to the physics of 
galaxy formation, and forms a readily observable quantity that can 
be compared with simulations. 

Numerous physical processes affect the evolution of a galaxy 
as a function of its density environment. The debate is ongoing as to 
which of these are intrinsic to the galaxies at formation and which 
are the result of later evolution. The relative importance of different 
processes depends on the density of the environment. For example, 
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in clusters infalling galaxies are subject to ram-pressure stripping 
and galaxy harassment. In groups, mergers, interactions and stran- 
gulation are more important. It is thought that ellipticals formed 
their stars early, exhausting their gas supply and are now passively 
evolving in dense environments, whereas spirals formed their stars 
more slowly, continuing to the present day and becoming subject 
to the physical processes described on infall into denser environ- 
ments as the structure evolved. We do not attempt to provide a 
full reference list for physical processes here, but rec e nt rev iews 
of the subject inclu de iBoselH & Gavazzil (l2006h . lHog j (l2006h and 
lAvila-Reesel ( 120071) . 

However, the morphology of the galaxies is affected in ways 
not necessarily correlated with the effect on spectral properties or 
colour. One can therefore investigate the morphology-density and 
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colour-density relations and subtract the two to see if a residual 
relation remains. If so, then there is an intrinsic effect of density 
that is not a result of a more fundamental process which causes the 
other relation. 

A particular example of this is the hypothesis that the proper- 
ties of galaxies within a dark matter halo depend only on the mass 
of the halo fe.g. lCoorav & S heth 2002) This can be tested both for 
colour and for morphology, although here we do not specifically 
relat e our results to haloes. Many results support this assumption 
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5007). 



Numerous existing results support colour as being a more 
fundamental predictor of environnient than morphology (e.^ 
KaufFmann et aDl2004l: liTanton et al.ll20 05bl; iQ uintero et al.ll2005 



Martinez & Muriell2006l : IOuintero et al . 2006). Ho wever , the ques- 



tion is not yet necessarily settled as Park et al. ('2007') find that 
the morphology and luminosity are more fundamental. Similarly, 
Lane et all feOOVl) study the Abell 901/902 supercluster and find 
that the local en vironment is niore im portant than the global for 
morphology, but lEinasto et "Zl (l2007h st udy numerous super clus- 
ters in the 2dF Galaxy Redshift survey ( Colless et al. "2001") and 
find that both the cluster and supercluster environment influence 
morphology. 

An example of a physical proc ess causing a distinct change in 
morphology is that of SO galaxies. IChristlein & Zabludofj (l2004h 
find that they must hav e formed from bulge enhancement and not 
disc fading. Bosell i & Gavazzil (f2006) find that SOs are formed 
through gravitational int eraction rather than spiral interaction with 
the intergalactic medium. [Moran et aD ( l2007h directly image spirals 
in the process of being transformed into SOs. 

In the S PSS, the colour-den sity relation was previously in- 
vestigated by iBalogh et aP (12004 hereafter B04). They found a 
pronounced bimodal distribution well-fit by a sum of two Gaus- 
sians. At fixed luminosity, the mean colours of the two distribu- 
tions were approximately independent of environment, but the red 
fraction strongly increased with density. The density measure used 
was the projected distance to the fifth nearest neighbour in a red- 
shift slice of width 1000 km . The part of our work dealing with 
u — r significantly parallels this work, although for a sample size 
containing three times as many galaxies. 

Previou s work in the SP SS on the morphology-density rela- 
tion includes Rjoto et al.l(l2003h . in which they use the conc entration 
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index and the texture parameter of lYamauchi et al. 
galax i es in the SPSS Early Pata Release (EPR, 
l2002h . ISorrentino et all ([2OO6) use the - r colour, the eClass 
eigenclass spectral type and the fraction of the galaxy profile de- 
scribed by a de Vaucouleurs profile (FracPeV) to describe their 
sample, with an en iphasis on comparing galaxies in voids to the 
general population. IPark et al.l ((2007) use an adaptive kernel esti- 
mate of density and the mo rphological separation via colour gra- 
dient of IPark & Choil (l2005l) and find little residual dependence of 
properties on environment at fixed luminosity and morphology, ex- 
cept for a sharp decrease in late types at the bright end. Some of 
these results are described further in ^3. 5 1 below. 

Here, we use the Fourth Pata Release (PR4, 
lAdelman-McCarthv et al. | l2006h of the Sloan Pigital Sky Survey 
fSPSS. Ivork et al.l 120001) to investigate the colour-density and 



morphology-density relations in a large volume-limited sample. 
We divide the sample by density and luminosity and investigate 
the trends of colour and morphology as a function of each, fitting 
Gaussians to the bimodal distributions in colour and morphology. 
The sample size, quality of photometry and corresponding 



spectra enable robust conclusions to be drawn from the data. How- 
ever, the available resolution limits us to simple measures of colour 
and morphology. Accordingly, as done by previous authors, we use 
the 1^ — r colour, inverse concentration index, C/inv and the Sersic 
index, n, each single-number axisymmetric measures of galaxy 
properties. However, in addition to this, we supplement the mor- 
phological measures w ith the Hubble type assigned by an artificial 
neural network fANN. lBall et al.ll200 4). This assigns T types 0-6, 
corresponding to E, SO, Sa, Sb, Sc, Sd and Im, to a limiting mag- 
nitude of r < 15.9 (compared to r < 17.77 for — r, C/inv and 
Sersic n). 

The Hubble type utilises more information from an image, for 
example spiral arm structure, and as such may show patterns un- 
seen in the C/inv and Sersic n. G iven the wide usage over time of 
this measure fe.g. lSandagell2005l) it is important that its usefulness 
is investigated in this context: if the results using the simple Hubble 
type used here are not useful even with the quality and quantity of 
the data we utilise, this will emphasize the requirement for differ- 
ent and more detailed measures of morphology for progress to be 
made. 

This sample used here is the largest to date for which the Hub- 
ble type measure of mor phology has bee n investigated as a func- 
tion of environment (c.f. Ivan den Berghl 120071) . who use the de- 
tailed eyeball- as signed Hubble type s for 1,246 galaxies from th e 
Rev ised Shapley-Am es Catalogue iSandage & Tammam] Il98ll) : 
and Conselice (2006), who construct a sample of 22,121 galax- 
ies from the Third Refere nce Catalog of Bright Galaxies (RC3, 
Ide Vaucouleurs et al.lll99lh . including Hubble types, but are inves- 
tigating ensemble properties rather than environment, which would 
require a more restricted sample due in part to the inhomogeneity 
of the RC3). 

Throughout, the standard spatial geometry is assumed, with 
Euclidean space, llmatter — 0.3, Ha = 0.7 and dimensionless 
Hubble constant h — 1, where h — ifo/100 km Mpc~^. 



2 DATA 

The SPSS is a project to map tt steradians of the northern galac- 
tic cap in five bands {u, g, r, i and z) from 3,500-8,900 A. This 
will provide photometry for of order 5 x 10^ galaxies. A multi- 
fibre spectrograph will provide redshifts and spectra for approxi- 
mately 10^ of these. A technical summary of the survey is given in 
lYorketal.1 12000). The telescope is d escribed in.Gunn et al. ( 2006) . 
The imaging camera is described in iGunn et al.l ( 1998 | ). The pho- 
tometric system and calibratio n are describ e d in [ Fukugita et al.| 
( 1996), Hogg et al. (2001), Smit h et al.1 (l2002h . llvezic et al. (2004) 
and Tucker et al .1(1200 6). The astrometric calibratio n is inlPier et al. 
(2003) and the data pipelines are in Lupton etliD (I2OO1I) . iLuptonI 
(2007) for the deblender, Frieman et al. and Schlegel et al. (in 

preparation). 

The targeting pipeline JStrauss "^^^[2002!) choo ses targets for 
spect roscopy from the imaging. A tiling algorithm dBlanton et al.l 
l2003h then assigns the spectroscopic fibres to the targets, the main 
source of incompleteness being the minimum distance of 55 arcsec 
between the fibres. This causes about 6% of galaxies to be missed; 
those that are will be biased towards regions with a high surface 
density of galaxies. The algorithm gives a more uniform complete- 
ness on the sky than a uniform tiling by taking into account large 
scale structure, but the effect is still present. The other main source 
of incompleteness is galaxies blended with saturated stars, which is 
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a 1% level effect. The overall s pectroscopic comp leteness is there- 
fore estimated to be over 90% (IStrauss et alJl2002h . 

The SDSS galaxies with spectra consist of a flux-l imited sam- 
ple, Main, with a median redshift of 0.104 dStrauss e t al. 2002); 
a Luminous Red G alaxy sample (LRG), a pproximately volume- 
limited to z ^ 0.4 (lEisenstein et al.ir200lh : and a quasar sample 
teichards et al.ll2002h . The limiting magnitude for the Main spec- 
tra is r < 17.77, which is substantially brighter than that for the 
imaging so the redshift completeness is almost 100%. A typical 
signal-to-noise value is > 4 per pixel and the spectral resolution is 
1800. The redshifts have an RMS accuracy of ±30 km s~\ 

We use galaxies from the Main Galaxy Sample. The data 
are extr acted from the New York Value-Added Galaxy Catalogue 
(VAGC, iBlanton et al.]r2005 a). They are maske d, extinction- and 
K-corrected in the same way as the datasets in iBall et all fcOOd 
hereafter B06). There is no correction for evolution, but as the red- 
shifts are limited to 2; < 0.1 its effects are not large. Similarly, 
there is no correction for dust, but this is mitigated by the require- 
ment that the galaxy axis ratio be less than that of an E7 galaxy, as 
described in B06. 



2.1 Galaxy Environment 

The masking of the dataset as in B06 leaves a sam- 
ple of 489,123 galaxies. These are matched to the galax- 
ies in the Pittsburgh-CMU Value Added Catalogue (VAC, 
http://nvogre.phyast.pitt.edu/dr4_value_added). The matches were 
carried out using the unique identification, both in object and date 
of observation, provided by the spectroscopic plate, mjd and fi- 
bre values for each object. The matches were made in order to 
utilise their measures of the distance from each galaxy to the near- 
est survey edge, which use a sophisticated algorithm taking into 
account the full survey mask. Galaxies measured by the VAC as 
being nearer to a survey edge than their Nth nearest neighbour are 
excluded to prevent the densities being biased downward near the 
survey edge or at the higher redshifts in the sample. 

We do, however, calculate our own environmental densities, 
as the VAC density sample is restricted to the range 0.053 < z < 
0.093 to gain an improved measure of the star formation rate, which 
is not considered here. We follow the method of B04: the density 
is given by Eat = N/Tid%, where c/at is the distance to the A^th 
nearest neighbour within ±1000 km in redshift. This is used 
to minimize contamination from interlopers. Here the value of N 
used is 5, fol lowing B04 who choose this value to approximate 
iDresslej {1980) who uses N — 10 before correction for super- 
imposed galaxies. We bin the sample in steps of 0.5 over the range 
-1.5 < log Us < 1.5, a factor of 1,000 in density. 

We set a limiting absolute magnitude for galaxies to be 
counted as neighbours of Mr — 51og/i < —19.5, corresponding 
to one magnitude fainter than the value of M* found for the overall 
Schechter fit to the galaxy luminosity function of B06. The cor- 
responding limiting redshift, at which a galaxy of this magnitude 
has an apparent magnitude r = 17.77 is z < 0.0889. This gives a 
sample volume-limited to a distance comparable to the z < 0.08 of 
B04. We also require z > 0.001 to exclude very local objects. We 
measure densities for galaxies in the range M*+^ i.e. -22.5 < 
Mr — 5\ogh < —17.5, binning in steps of one magnitude. This 
again is similar to the range of -22.23 < Mr - 51og/i < -17.23 
of B04. 



2.2 Galaxy Colour and Morphology 

We use u — r colour following B04 and lBaldrv et al.l (l2004bl) . who 
find that the galaxy population is well described as the sum of 
two Gaussians i n u — r, with an optimal colour separator of 2.22 
( Strateva et al. 200ll). We use the model colours from the same 
dataset as B06, extinction- and K-corrected, as in that paper, to a 
rest-frame band- shift of z = 0.1, denoted throughout in the plots 
by {u — r)o.i, and in the text by — r. The model colours are de- 
rived from the SDSS model magnitudes. For each galaxy these are 
derived from the best fitting of a de Vaucou leurs Cde VaucouleursI 
Il948h or exponential profile (lFreemanlll97Qh . 

For morphology, we use three measures: the inverse concen- 
tration index, C/inv, the Sersic index, n and the Hubble type. The 
C/inv and Sersic index provide a measure of morphology that, 
while simple, is available for the full r < 17.77 sample. C/inv is 
directly from the SDSS DR4 and is given by R50/R90 where Rn 
is the NVo light Petrosian radius (see, e.g. IStoughton et al."2002S. 
The inverse is used because it has the range 0-1. The Sersic in- 
dex, from the VAGC, is the expo nent in the for mula for a gener- 
alised gala xy light profile given bv lSersid(ll968h . or more recently 
in e.g. Gra ham & Drivej (l2005l) . Because n is an exponent, and 
the bimodality shows more clearly when it is done so, n is binned 
logarithmically, as with the density. As with the colours, these are 
binned appropriately, giving bins of 0.23 < C/inv < 0.58 in steps 
of 0.07 and -0.25 < log n < 0.75 in steps of 0.2. We truncate 
the Sersic index at log n = 0.75 to avoid a spike in the number 
of galaxies in the VAGC data with log n ^ 0.77. Their data is 
also truncated at n = 6 (log n = 0.778). Parameters other than 
colours are measured in the r band, since this band is used to de- 
fine the aperture through which Petrosian flux is measured for all 
five bands. 

The ANN morp hologies u s ed we re also the same as those in 
B06, assigned as in lBalletal.1 (l2004l) and updated as described 
in B06. The types assigned correspond to the Hubble types E=0, 
S0=1, Sa=2, Sb=3, Sc=4, Sd=5 and Im=6. As described in B06, 
there is bias away from the ends of the scale, particularly for late 
types. This is probably due to the dominant contribution of the con- 
centration index in the training set, which shows a similar bias, and 
the relative lack of types of Sd or later in the training set, about 5% 
of the total ( Fukugit a et al.ll200 7) in part due to the r band selec- 
tion of the SDSS instead of the more commonly used bluer bands. 
Hence the latest types assigned have Tann ^ 4.5. Fig. [T] shows a 
random sample of five galaxies from each T type, rounded to the in- 
tegers 0-4. We bin the ANN morphologies in five bins from —0.5- 
4.5. The colours are also binned in five bins. Unlike the other mea- 
sures presented, the limiting magnitude for the ANN morphologies 
is r < 15.9, as this was the faintest level to which the visually 
assigned training set was generated. Of the galaxies studied here, 
13,655 have a value for Tann and a valid density. Combined with 
the absolute magnitude limit of Mr — 5\ogh < —17.5, these limits 
ensure that there is no extrapolation from the ANN training set in 
apparent or absolute magnitude. 

For each set of measurements, the galaxies are required to 
have values within the ranges of the bins of all quantities that are 
binned, i.e. density, luminosity, and one or both of colour and mor- 
phology. This ensures that each related figure is showing the same 
galaxies. The resulting samples consist of 79,553 galaxies for i^ — r, 
79,606 for C/inv, 75,861 for Sersic index and 13,655 for Tann- 
For the residuals in u — r and morphology, requiring both colour 
and morphology to be within the bin ranges reduces the sample 
sizes to 79,495, 75,753 and 13,638 respectively. 
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Figure 1. SDSS cutout images showing five random galaxies of each of the T-types 0-4 (rounded to the nearest integer) on respective rows. The T types 
correspond to Hubble types E, SO, Sa, Sb and Sc. This shows the typical quality of the data used in determining the SDSS object parameters used in training 
the artificial neural network to assign the types. A larger selection of similar images can be found in lFukugita et al.l ( I2007n . 



2.3 Data Fitting 

The data points for the histograms of colour and morphology for 
bins of log S5 and Mr are fitted with a double Gaussian using a 
simplex search algorithm over the ranges 0.5 < u — r < 4.0, 
0.23 < C/inv < 0.58, -0.25 < n< 0.75 and < Tann < 4.5. 
The number of galaxies outside these ranges is small, being less 
than 1% of the total. The fits involve six parameters and are given 
by 



(1) 



where the N(x,a) are Gaussians with mean x and standard devia- 
tion a and A and B are scalars. Whereas B04 constrain the a values 



of the Gaussians to be functions of Mr only, here we do not apply 
any restrictions on the parameters of the two Gaussians in one plot 
relative to any others. The fits do not utilise the extreme ends of the 
distributions, where the number of galaxies is negligible. 



Again following B04, the la error bars for each bin are given 
by V A^bin + 2 where A^bin is the number of galaxies in that bin. 
This is an approximation of the Poisson error for small number 
statistics. 
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Figure 2. Contour plot of Petrosian inverse concentration index, Clmv, versus rest-frame u — r colour, subdivided into panels by environmental density, E5 
and absolute magnitude, Mr. The sample consists of 79,606 galaxies. The contours are normalised to the numbers of galaxies in each panel, so small numbers 
lie outside the lowest contour. 



3 RESULTS 

3.1 Two-Dimensional Distributions 

Figs. [JHH show colour versus morphology (C/inv, Sersic n and 
Tann respectively) for the sample, divided into panels according to 
density and luminosity. It is clear that there is a broad division into 
two populations. This division justifies the treatment below of the 
galaxies as two populations modelled by a double Gaussian, which 
is used in ^3.2j - ^3.4l below to quantify the trends seen in more de- 
tail. 

In Fig. [21 the populations are approximately centred on — 
r - 2.75, C/inv - 0.35 andu-r ^ 1.5, C/inv - 0.45. These 
correspond to the red, early-type and blue, late-types. As expected, 
there are few red galaxies with high C/inv or blue galaxies with 
low C/inv, and there is little dependence of each population centre 
on luminosity or density. What does change radically is the rela- 
tive numbers of galaxies in each population, which is emphasized 
here by the normalisation applied to the contours in each panel. 
Blue colours and diffuse morphologies dominate in galaxies of low 
luminosity in low density environments. They are almost entirely 
absent in luminous galaxies and in high-density regions. 



Fig. [3] shows a similar pattern, although the sense of n is in- 
verted compared to C/inv The truncation in the Sersic index at 
log n < 0.75 f ^2.2l ) is visible, although the contours show that few 
galaxies are excluded. At fainter magnitudes, Mr ^ —19.5, the 
visual overlap between the populations is lower in n than in C/inv 

Fig.|4]shows that there are again similar patterns in Tann mor- 
phology. Types < Tann ^ 2, i.e., E-Sa, are all approximately 
the same colour, although the incidence of types Tann ^ 1 de- 
creases at bright magnitudes Mr ^ —20.5. There are relatively few 
galaxies at the faint end, —17.5 < Mr < —18.5 due to the brighter 
r < 15.9 magnitude limit for Tann compared to r < 17.77 for 
C/inv and n. 



3.2 Gaussian Fits 

Fig.[5]shows the rest-frame u — r colour distribution, subdivided by 
density log S5 and luminosity Mr. The distributions are clearly bi- 
modal and well fit by a sum of two Gaussians, confirming the idea 
of two separate pop ulations described by numerous authors (e.g. 
iBaldrv et aDl2004bh . As expected, the fraction of blue galaxies de- 
creases with density. The position of the blue peak moves redwards 
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Figure 3. As Figure |2] but with C/inv replaced by the logarithmic Sersic index, log n, for 75,861 galaxies. The truncation of the data at log n = 0.75 ( ^22\ 
is visible. 



from u — r ^ 1.5toi^ — r ^ 2.25 and that of the red peak stays ap- 
proximately constant, with just a slight reddening from u — r ^ 2.5 
to 1^ — r ~ 2.75. Unlike the peaks, the position of the minimum 
between the two distributions stays constant, although in the higher 
density and luminosity bins the exact position of the minimum is 
unclear. Versus luminosity, a similar trend is seen with the red peak 
moving slightly bluewards and dropping in fraction with decreas- 
ing luminosity and the blue peak becoming bluer. Bright galaxies 
at Mr — 51og/i < —21.5 are relatively rare as are faint galaxies due 
to the r < 17.77 limit for spectra. 

The figure represents an update of fig. 1 of B04, using DR4 
rather than DRl and with the differences described in 32]above. As 
mentioned, B04 find that at fixed luminosity the mean colours of 
the two distributions are nearly independent of density, with just 
a weak trend for colours to become redder at high densities, but 
that the fraction of red galaxies rises sharply at all luminosities. As 
expected, our plot shows the same trends. 

Figs. [6H1] show the same as Fig. [5] but for the C/inv, 
Sersic n and Tann morphologies respectively, i.e., they show the 
morphology-density relation subdivided by luminosity. 

As with colour, the concentration index shows clearly bimodal 
distribution, although significant overlap between the populations 



causes a unimodal appearance in some of the panels. The plots 
show peaks at C/inv ^ 0.3 and C/inv ^ 0.43, corresponding 
to the de Vaucouleurs elliptical galaxy profile and the exponential 
spiral disc profiles, respectively. As expected, early-type galaxies 
dominate at high density and high luminosity. 

The Sersic index shows a bimodal distribution at all densities 
and all but the highest and lowest luminosities. The indices n — A 
(log n = 0.60) and n = 1 (log n = 0) correspond to the de 
Vaucouleurs and exponential profiles. The peaks in the distribution 
are close to these values, but slightly nearer the overall mean, in a 
similar manner to C/inv 

For Tann, a sum of two Gaussians is now a worse fit to the 
data, with a clear excess of galaxies at Tann 1 at densities 



log 



> 



-0.3. The expected dominance of early types at high 



density is seen and the peak narrows with density. The late-type 
peaks becomes lower and earlier, in a similar manner to the red- 
dening of the blue peak in — r. The lack of (or less clear) division 
into a bimodal population suggests that if galaxies really are di- 
vided into just two dominant populations then either (1) the Tann 
is insensitive to it, spreading between the two; (2) that Tann is 
more sensitive to transitional states between the two, such as SO 
galaxies, than colour, or (3) the Tann is an intrinsically 'fuzzy' 
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measure, resulting in a inevitable spread across types. If (2) is cor- 
rect, such a transitional state must then be rare and/or brief in colour 
space for the clear division in colour to be seen as SO galaxies are 
not rare. That the transitional population is clearest at higher den- 
sities is also consistent with it being SOs. (3) is less likely for the 
overall population, i.e. one can split E/SO from spiral, but the ex- 
cess over Gaussian at Tann ~ 1 and at very late types supports 
the 'fuzzy' nature of the Hubble type classification. The spread is 
not likely to be a n artefact of typing v ia neural network, as these 
were shown (e. g. TStorrie-Lombardi et al. 1992; Naimetal. 1995; 
iBall et al.ll200 4) to perform with an RMS spread between true type 
and assigned type no different to that between human classifiers. 



In most of the plots in Figs. [SHSl the reduced values, ex- 
pected to be around 1 or less, of the double Gaussian fit are large. 
This was also seen in the for fitting the bivariate LF in B06. As 
with those results, this is probably due to the small size of the error 
bars resulting from the large number of galaxies in the sample, indi- 
cating that the fitting function used is only an approximation to the 
data. What is clear is that the data is better fitted by two Gaussians 
than, say, one. 



3.3 Red Fraction, Mean, Dispersion 

Figs. l9tiTT] show some of the trends in Figs. [5HS1 more clearly. The 
±lcr error bars are calculated by finding the point at which the 
for the fit increases by 1 around the minimum for each parameter 
in turn. The shape of the is assumed to be parabolic, which is a 
good approximation, although it assumes that the is minimized 
for the other parameters whilst the one under study is being varied. 
However, these approximations will tend to overestimate the errors, 
so it is good as a conservative estimate. 

The top left panel of Fig. |9] shows the fraction of red galaxies 
for Mr versus for — r. The fraction is given by the ratio of 
the areas under the two Gaussians, which is given by A/{A-\- B) 
assuming that approximately the whole area under each Gaussian 
is populated by galaxies (the normalised areas are 1). This is seen 
to be true in Figs.|5}|8] The fraction rises with density, as expected, 
and rises in a similar way over all luminosities, although the bright- 
est galaxies at —22.5 < Mr — 51og/i < —21.5 are slightly noisy. 
The trends are similar to those in the left hand panel of fig. 2 in 
B04. 

The top left two panels of Fig. [Tol show the mean u — r colours 
xi and X2 in a similar way for the two Gaussians. The mean colour 
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Figure 5. Numbers of galaxies for w — r subdivided by E5 and Mr. The points are the number of galaxies in each bin; the error bars are given by V A/'bin + 2 
where A/'bin is the number of galaxies in that bin; the solid line is the best fit of the sum of two Gaussians; the dot-dash and dashed lines are the individual 
Gaussians. The vertical axes within each row of panels are given equal heights for clarity. This plot is similar to that of .Balogh et al.. (,20Q4i) . but for DR4 as 
opposed to DRl. The trends seen are very similar, as expected, because DR4 is a superset of DRl. 



for the red galaxies is approximately independent of density (as 
seen in Fig.|5]), with just a slight trend for redder colours at higher 
densities. The blue galaxies show a more distinct reddening to- 
wards the highest densities. This is similar to fig. 3 of B04. 

The top left hand two panels of Fig. [TT] show the — r disper- 
sions cTi and (72. For red galaxies the dispersions are low and for 
lower luminosities decline at the highest density, whilst those for 
the blue galaxies are higher, slightly increasing at the highest den- 
sity for low luminosity. The dispersions for red galaxies also show 
a clear decrease with increasing luminosity. 

The remaining panels of Figs. I9HTT] show the same as — r 
but for the C/inv, Sersic n and Tann morphologies. Each measure 
of morphology shows a similar behaviour. In Fig.[9j the early-type 
fraction shows some increase with density in a similar way to the 
red fraction, whereas the late-type fractions (not shown) are noisy. 

In Fig. [To] the mean type for the early-type galaxies is seen 
to be approximately independent of density at fixed luminosity in 
the same way that colour was for the red galaxies. This hints that 
the same population is being seen, which is consistent with a dis- 
tinct bright, concentrated, early spectral type, red population. The 



late types change little with density, the C/inv and Tann becom- 
ing earlier but the Sersic n showing no trend. The similarity of the 
patterns seen suggests that in terms of the effect of environmental 
density, colour is closely correlated with morphology, but that there 
could be a residual relation. 

Fig.[TT]shows that for morphology the dispersions do not show 
any obvious trends with density. For u — r, this justifies the con- 
straint in B04 for the the dispersion not to vary. (B04 also find little 
change if they relax this constraint.) Some of the larger changes are 
due to the widths of individual Gaussians being poorly constrained 
by either significant overlap between them or by low numbers of 
galaxies. 



3.4 Residual Colour- and Morphology-Density 

The residual relations are calculated by binning the galaxies into 
a 2-dimensional grid of colour and morphology, using the same 
histogram bins as in Figs.|5}|8] The mean density in each bin is 
calculated, along with the marginal averages as a function of colour 
bin (/ij) and of morphology bin (/ij^), using the densities for the in- 
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Figure 6. As Fig.[5]but for CIu 



dividual galaxies in each so that the bins are correctly weighted. In 
order to factor out the colour-dependence of density, one divides 
each mean density /j^ij by the appropriate marginal average /i^. 
A re-calculation of the marginal dependence of density on colour 
would then give all ones, whereas a re-calculation of the marginal 
dependence on morphology will yield the residual morphology- 
dependence of density after that due to colour has been removed. 
By instead dividing out the morphology-dependence in an equiv- 
alent way, the residual colour dependence may be calculated. The 
errors on each point are the la dispersion of the values of the bins 
which are averaged to give that point. 

Fig. [12] shows the residual relation for — r on removal of 
morphology. There is a definite trend in density with colour, con- 
sistent with numerou s earlier studies (e.g. Kauffma nn et "^1 120041 : 
iBlanton et alj|2005bh in which the colour was a strong predictor of 
environmental density. As expected, the density increases for red- 
der galaxies. The error bars are still relatively large, reflecting the 
less certain nature of S5 compared to colour. 

Fig. [13] shows the morphology-density relation after the re- 
moval of the change in density due to the change in colour with 
morphology. In contrast to colour, the residual relation is consis- 
tent with a value of one, i.e. no residual relation. This means that 



there is no evidence for a change in morphology that is not corre- 
lated with a change in colour. 

3.5 Comparison to Previous Work 

Besides B04, who presented the colour-density relation for the 
EDR and for whose results ours broadly represent an update to the 
SDSS DR4 ( ^331 and ^Mlabove), substantial work involving simi- 
lar measures of colour and morpho l ogy w ith comparison to density 
has been carried out bylGoto et alJ ([2003') and Park et all i2007h . 

iGoto etaD(l2003l) present the morphology-density relation for 
the SDSS EDR. They assign morphology using the inverse concen- 
tration index, as used here, and a preliminary version of the texture 
parameter described by Yamauchi et al. (2005). The latter is cor- 
rected for elliptical isophotes and gives an improved correlation 
to visual morphology. The galaxies are divided into early, inter- 
mediate, early disc and late disc. The sample is volume-limited to 
Mr — 51og/i < —19.9 in the redshift range 0.05 < z < 1 and con- 
tains 7,938 galaxies. This compares to our values of Mr — 51og/i < 
-19.5, 0.001 < z < 0.0889 and 13,655 galaxies. The same 
^matter = 0.3, Q A = 0.7 Euclidcan cosmology is used. Their 
density estimator is three dimensional. 

They find the expected trend of increasing early-type fraction 
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Figure 7. As Fig.|5]but for Sersic \ 



with increasing density, and also two characteristic scales, giving 
three density regimes. In the lowest density regime, the relation 
is less noticeable. In the intermediate density regime, intermedi- 
ate type galaxies increase with density and spirals decrease. In the 
densest regime, the intermediate types decrease and the early types 
increase. Thus the morphology is little affected in the low density 
regime, then subject to a mechanism which inhibits star formation 
and turns spirals into intermediate types in the intermediate density 
regime, then a third process decreases these types and increases 
early types in the dense regime. Our work is consistent with this, 
with the bridge between the two bimodal populations resembling 
the intermediate types described here. It is difficult to compare 
quantitatively, however, as the density measures are different. 

|ParketaD(l2007h study numerous galaxy properties as a func- 
tion of environment in various volume-limited samples containing 
up to 80,479 galaxies for 0.025 < z < 0.11 using a subset of the 
SDSS DR5. The galaxy properties studied include the {u — r)o.i 
co lour an d morphology assigned using the colour gradient method 
of IPark & Choi (2005). Their measure of density is obtained using 
an adaptive smoothing kernel, so again an exact quantitative com- 
parison with our results would be difficult. Their fig. 9 shows th e 
fraction of early-type galaxies (as defined in IPark& Choil ( l2005h ) 
as a function of density and luminosity. The result shows the same 



trends as our Fig. 5, with a monotonic increase in early-type frac- 
tion with density and luminosity, and the trend lines separated by 
luminosity approximately parallel to each other. However, they also 
fix both luminosity and morphology and then find little residual de- 
pendence of properties on environment, except for a sharp decrease 
in late types at the bright end. 



4 DISCUSSION 

Considering the number of galaxies, the variety of their properties 
and environments and the quality of the SDSS data, the rest-frame 
colour histograms are well fit by a simple sum of two Gaussians. 
This result is also seen by B04 and others and strongly suggests 
two underlying populations of galaxies. This is consistent with the 
patterns seen in the LPs in B06 (their figs. 2, 3, 4 and 6 show the LF 
bivariate with Tann, C/inv, Sersic n and (u — r)o.i respectively). 
There, an early morphological type, concentrated, early spectral 
type, red, high Sersic index population of galaxies is evident and 
the LPs binned in absolute magnitude and plotted against their sec- 
ond parameter often show bimodality. 

One of the main questions asked in this paper is whether the 
morphology reveals anything that the colours do not. Compared to 
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Figure 8. As Fig.|5]but for Tann- 



u — r, the C/inv and Sersic n also show clearly bimodal distribu- 
tions, although the C/inv plot becomes unimodal at fainter mag- 
nitudes and low densities. There are no obvious residual relations 
when those due to colour are removed, thus these measures of mor- 
phology appear not to be predictive of environment beyond their 
correlation to colour. 

When the Tann morphology-density relation is plotted the 
distributions are less bimodal, ranging from looking similar to the 
plots for colour to being not bimodal at all. There are various possi- 
ble reasons for this, including (1) the morphology may be less sen- 
sitive to the processes producing the bimodal population in colour, 
being intrinsically spread when compared to the colour-producing 
physics; (2) there may be more than two populations, the third or 
higher to which colour may not be sensitive (an example might be 
SO galaxies, which are almost as red as ellipticals); or (3) the neu- 
ral network morphological types might be spread or bia sed versus 
the tr ue types. The last point is not thought to be likely jBall et al.l 
I2OO 4Y It is difficult to distinguish between the intrinsic spread or 
more populations possibilities, although the fact that some of the 
bins are more bimodal than others favours the third population idea. 
Here neither possibility is ruled out and indeed both could be occur- 
ring. More detailed investigation of the morphology-colour-density 
relationship may be useful, but it may remain the case that the mor- 



phologies are simply intrinsically too 'fuzzy' to show the underly- 
ing populations as clearly as the colour. That said, Tann has larger 
error bars, so the formal goodness of fit is not necessarily worse 
than that for — r, C/inv and n. 

This spreading of Hubble types between the two populations is 
also seen by vari ous o ther worker s. In the Millennium Galaxy Cat- 
alogue fMGC. lLiske et alJl2003h JPriver et all feood) study a sam- 
ple of 3314 galaxies with eyeball Hubble types at B < 19 and find 
that, while types E/SO (bulges) exclusively occupy the early-type 
peak and Sd and later (discs) the late-type peak, the Sa-Sc, which 
are a mixture of b ulge and disc, are spread between the two peaks. 
Also in the MGC jEUis et all (l2005h study a well-defined sample of 
350 galaxies and find, using multivariate statistical analysis of the 
galaxy properties, that only two populations, corresponding to early 
and late types, are justified, and the late types (Sa-Irr) are smoothly 
spread within the late-type region of parameter space. In a sample 
of 1,246 galaxie s with Hubble types from the Revised Shapely- 
Ames catalogue, van den Bergh" uOoT) similarly find a continuum 
of Hubble types but a bimodality in colour. 

In a sample of 22,121 gala xies from th e Third Reference Cata- 
logue of Bright Galaxies, .Conselicd (I2OO6I) finds that the best three 
parameters with which to classify galaxies are mass, star formation 
and merger history. The Hubble type correlates most strongly with 
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Figure 9. Red and early-type galaxy fractions in each Mr bin versus E5 for w — r, C/inv, Sersic n and Tann- The error bars are Icr. The left hand panel 
shows a similar upward trend to fig. 2 of B 04. 



stellar mass, {U — B)o and {B — V)o colour. They also (his fig. 
10) confirm the intrinsic spread in morphology when compared to 
other properties such as colour, showing that morphology is at least 
to some extent independent of the other parameters. 

iBaldrv et al.l (l2004ab discuss the work in lBaldrv et aP (l2004bh 
and B04 and agree with the interpretation of morphology, in the 
sense that 'whatever processes give rise to the blue/red distribution 
should also give rise to distributions in morphology', although they 
go on to say that 'thus, SO or Sa galaxies could have a probability 
of belonging to one or the other distributions and should not be 
considered as classes'. 

The physics giving the bimodal population in colour is thought 
to be that red galaxies have a passively evolving old stellar popu- 
lation and are formed by major mergers. The blue galaxies are un- 
dergoing star formation and less violent accretion events. Galaxies 
can change from blue to red and when they do they do so rapidly, 
leaving few galaxies in between the two populations. B04 go on to 
propose that 'most star forming galaxies today evolve at a rate that 
is determined primarily by their intrinsic properties, and indepen- 
dent of environment', which is consistent with the results here. The 
morphology may be seeing extra physics at work, for example ram- 
pressure stripping of gas from a galaxy infalling into a cluster can 
turn a spiral into an SO rather than an elliptical, even though they are 
almost as red. SOs are sufficiently common that they should show 



up as a component in the plots in this paper if they indeed form 
a component distinct from other galaxies. In comparison to simu- 
lations, only recently have semi-analytic models be gun to produce 
bimodal populations matching those observed (e.g. ICattaneo et al.l 
12006^, so there may still be much to be gained by a detailed inter- 
comparison between the two. 

The ANN types are biased away from very early or late types, 
but not to an extent large enough to affect the overall trends seen 
here. This is thought to be due to the importance of the concentra- 
tion index in the training set, which has a similar bias. A larger 
training set for morphology is needed, particularly for very late 
types of type Scd or later, corresponding to Tann ^ 4.5. Our 
training set was based on the SDSS EDR version of the catalogue 
of Fukug ita e t al. ( 2007). However, their published version, while 
based on DR3 and 40% larger at 2658 galaxies instead of 1875, 
would be unlikely to significantly improve our results due to the 
intrinsic spread in Tann and the r-band object selection of the 
SDSS, which limits the Scd or later types to 5.2% of the sam- 
ple. They compare th is to 10.6% from the earlier B-band study of 
lFukugitaeta"lT(ll998h . 

Here the model u — r colour was used. [Driver et al ] (l2006h 
showed that the core u — r, measured using the PSF colour, pro- 
vides a cleaner separation for their sample when used in conjunc- 
tion with the Sersic index. The same may be true here if that u — r 
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were used, although the model u — r is cl eaner than the Petr osian 
u — r, for which we also generated results. IPriver et al] (l2006h also 
argue that the two fundamental components are, rather than red and 
blue galaxies, bulges and discs. This is consistent with the results 
here. H owever, it is not necessarily that simple, as Drorv & Fisher 
toot differentiate classical and pseudobulges, which have differ- 
ent formation mechanisms, and find that, for types Sa-Sbc for a 
given bulge-to-total light ratio, the former are predominantly in red 
galaxies and the latter in blue, i.e., the bimodal distribution is not 
simply due to the bulge to disk ratio. 



5 CONCLUSIONS 

Galaxy properties are studied as a function of environment in the 
SDSS, specifically the relations between u — r restframe colour, 
morphological type and the environmental density. The colour is 
given by the restframe value of u—r. The morphology is in the form 
of the inverse concentration index, C/inv, the Sersic index, n and 
the Hubble T type, Tann, the latter assigned by a trained artificial 
neural network. The density is the Ss surface density of galaxies 
measured by the fifth nearest neighbour brighter than an r band 
absolute magnitude of Mr — 51og/i = —19.5. The range of density 
probed covers all environments from the cores of rich clusters to the 
field. The colour results are similar to those of lBalogh etaLl ( l2004h 
(B04) but here these are extended to galaxy morphology. 



Subdivision of the population by density and luminosity 
shows a clear bimodal distribution in the colour-morphology plane. 
For each subsample we fit the colour and morphology by a sum of 
two Gaussians. The well-known colour-density and morphology- 
density relations are seen. The colour is well fit by the Gaussian 
sum, confirming two underlying populations, one red and one blue, 
as seen by B04. The mean colour of the red population is approx- 
imately constant but the blue population gets slightly redder with 
increasing density. The relations are seen at high significance over 
the range of densities probed, thus confirming that they extend into 
the field and are not just present in galaxy clusters. 

Compared to — r, C/inv and n, the Tann morphology is 
less obviously fit by the sum of two Gaussians, suggesting one or 
more of the possibilities (1) three or more populations (e.g. early, 
lenticular and late), (2) sensitivity of the morphology to processes 
which do not affect the colours, or (3) a higher intrinsic spread in 
morphology compared to two underlying populations and physical 
processes. 

While the C/inv and Sersic n to some extent separate the two 
main galaxy populations, certainly more cleanly than the Hubble 
type, both still suffer from extensive overlap in the populations. 
Figs. [2H4] show that, rather than using a single colour or morphol- 
ogy, improved separations should be obtained using a division in 
the plane of colour and morphology rather than any single measure. 
In its simplest form in Figs. [JHJl this would correspond to the al- 
lowance of an oblique line as a separator instead of an axis-parallel 



© 2007 RAS, MNRAS 000,[TJtT6l 



14 A^. M. Ball et al 




line (c .f. e.g. fig. 6 of IChoi et al l (l2007h . fig. 10 of iBaldrv et aP 
([2006^, or fig. 15 of IPriver et all (l2006h . which also show this pla- 
nar distribution). 

The colour-density relation divided into bins by absolute mag- 
nitude reproduces the results of B04: the mean colours are inde- 
pendent of density for fixed luminosity but the red fraction in- 
creases with increasing density. The fraction of red galaxies in- 
creases approximately linearly with density, with parallel slopes of 
lower fraction for decreasing luminosity. The dispersion in colour 
for each population is low for red galaxies, increasing with decreas- 
ing luminosity and decreasing with density for low luminosity. For 
blue galaxies it is significantly higher increasing at low luminos- 
ity with density. At high luminosity there is now a significant trend 
with density. 

The morphology-density relations divided in this way are sim- 
ilar if the early-type galaxies are red. The early-type galaxy frac- 
tions increase in a similar way to the fraction of red galaxies but 
the trends are noisier. The mean types versus density are similar to 
colour but the late types do not as obviously get earlier with density, 
consistent with the colour having a stronger dependence on density 
than morphology. The dispersions show few significant trends with 
density. 

Inu — r there is a clear residual relation of galaxies for both 
blue and red galaxies with density even after the reddening due to 
the increase in early types is removed. This suggests that colour has 
a stronger dependence on density than morphology, in agreement 
with other studies. 

For each of C/inv, n and Tann, there is no evidence of a 



residual trend once that due to variation with colour is removed. 
Therefore, for these data, none of the measures of morphology are 
predictive of environment beyond their correlation to colour. 
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Figure 12. Residual density dependence of colour once the dependence due 
to morphology (C/inv. Sersic n and Tann in respective panels) has been 
removed. The sample sizes are 79,495, 75,753 and 13,638 galaxies respec- 
tively. The error bars are the la dispersion in the morphology bins, and a 
value of 1, highlighted by the horizontal line, indicates no residual. There is 
a clear reddening of colour with density. 
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